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For transport model simulations of collisions between two nuclei which have N/Z significantly 
different from unity one needs a one-body potential which is both isospin and momentum dependent. 
This work provides sets of such potentials. 
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I. INTRODUCTION 
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^ Momentum dependent mean-fields for transport model calculations of heavy ion collisions have been in usage for 

/— s ' many years now [1-6] . So far the attention has been for a momentum dependent potential which does not distinguish 
between neutrons and protons. This is adequate for systems which have N k, Z where N is the number of neutrons 
\ and Z the number of protons. One main focus of nuclear physics research today is to explore both structures of 
, exotic nuclei in regions far-off the stability line and novel properties of neutron-rich nuclear matter. The latter can be 
investigated by using collisions induced by neutron-rich nuclei at intermediate to high energies. To interpret critically 
, data from these collisions and to extract accurately properties of neutron-rich nuclear matter, advanced transport 
model calculations are necessary. In asymmetric nuclear matter, the one body potential seen by a proton is different 
from that seen by a neutron. This has been implemented in BUU (Boltzmann-Uehling-Uhlenbeck) calculations but 
&^ with a simplification that the potentials are taken to be momentum independent [7] . The present work aims to correct 
C^l this deficiency. That the momentum dependence will be different for neutrons and protons is of course well-known 
^""^ . and has been the subject of quite sophisticated many body calculations, see, e.g, ref. [8] for a recent review. We do not 
' aim to add anything fundamental in this regard. Our objective is to obtain a parametrized version which displays the 
, main characteristics of momentum dependence in asymmetric matter and is still usable in practical BUU calculations. 
' Major advances in this direction were already made: see articles by Bombaci [8] and Prakash et al. [9]. We add to 
"*7^ \ this. We will not only extend the simplest momentum dependent potential [1] to include isospin but also extend the 
Q ■ improved treatment [2-4] subsequently introduced to include isospin. Thus this is an extension of the work reported 
^ ' in [3,4]. 
^ ■ 
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^ An effective momentum dependent potential can be deduced from phenomenological interactions. We take the 

■ " " ' Gogny interaction [10] to obtain an idea of the momentum dependence. There are many reasons for this choice. It has 
been used in detailed fits for spectra in finite nuclei. It gives accepted values for binding energy, saturation density, 
compressibility and symmetry energy in nuclear matter. It has been verified already [2] that the interaction produces 
a reasonable parameterisation for the real part of the optical potential in nuclear matter as a function of incident 
energy. The Skyrme interaction has a wrong asymptotic behaviour as a function of energy. (This is amplified in 
[2].) Since we want to devise a momentum dependence which should hold for beam energy as high as 1 GeV/nucleon 
(this would allow investigation of symmetry energy at higher than normal nuclear density) we discard the Skyrme 
interaction. 

For the purpose of this work we will define nuclear matter to be an infinite system but without the restriction 
N = Z. Using the Gogny interaction, we deduce U{p, S,p,t), the one body potential a particle of momentum p and 
isospin T feels in cold nuclear matter with density p and asymmetry S = y-^- One then generalises to U in the 
case of heavy ion collisions. For BUU calculations U is the only quantity needed. But it is useful to also have an 
expression for V{p,5), the potential energy density in cold nuclear matter with a given density p and asymmetry 5. 
This allows one to deduce E/A as a function of p and S which is, of course, of importance. The expression for V{p, S) 
can also be generalised to the case of heavy ion collisions and can be used to check, for example, the accuracy of 
energy conservation in a BUU simulation. 



II. A MOMENTUM DEPENDENT POTENTIAL FROM A PHENOMENOLOGICAL INTERACTION 
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Fig.l: Single particle potential U with respect to k and total single particle energy (e) for neutron and proton for 
6=0.4- The solid line is the single particle potential for symmetric matter; here p = 0.16/m~^. 



Normally nuclear matter denotes an infinite nucleus with N = Z. Total potential energy in cold matter is deduced 
from 



Vt = ^ 51 {Pi,(^i,ri,P2,(T2,T2\v{r) (2.1) 

pi(TlTip2CT2T2 

(|pi,ai,Tl,P2,0-2,T2)- \p2,<J2,T2,Pl,ai,Tl)) 

where, 

=^{W + BP„- HPr - MP„Pr)ie-'""^^'" (2.2) 

i=l,2 

+^o(l + P.)p"(^^^)^(r-'l-f2) 

There are two finite range Gaussians and a density dependent zero-range force. The values of the parameters are 
given in [10]. The one body potential U{p,S,p,T) is obtained from 

U{p,S,p,t) = X]p'cr'r'(P''^'''''P '''^''■''') ~ |P CT, t)) plus rearrangement term which 
for nuclear matter is {3/2)taap"-^{l/4:)p'^{l - S'^). 

The momentum dependence in U comes entirely from the exchange term of the finite range part, i.e., from 
{p,p '\e~'^ 1^ \p ' tP). Except for the momentum dependent part, very simple expressions for U and Vr/^, the 
potential energy per particle are obtained for the Gogny potential. Thus U{p,b,p,T) = X + Y + Z where X arises 
from the direct term of the finite range interaction, Y arises from the to term (density dependent two body term) and 
Z from the exchange term of the finite range interactions. For a given r, these are: 

^ = H E ^'/V?(W- + B/2)A - p. ( E ^'^Vf (^? + M/2)A (2.3) 
Y = ltop"{p-Pr) + ltoap^-'\p\l - 6') (2.4) 
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Here r' ^ r, the isospin of the particle whose one body potential is being sought. 

Similarly, Vt /A=potential energy per particle, has contributions from the direct term of the finite range force, from 
the density dependent term and the exchange term of the finite range force. Denoting the first two by X' and Y' 
respectively, explicit expressions for these are: 
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Fig. 2 Esym{p) as a function of p. 

Wo do not write down the explicit expression for Z'. It is obtained from the expression of Z above after a further 
integration over p, sum over a, r and dividing the answer by 2. 

The one body potentials as a function of k for neutrons and protons in cold matter as predicted by the Gogny 
potential for S = and S = 0.4 arc shown in Fig. 1. They arc quite similar to other calculations of U{p,5^p,T) 
available in the literature. We compare, in particular, to the U{p,d,p,T) given in Fig. 4 in [11]. That figure is for 
p = .17 and obtained from Brueckner-Hartee-Fock calculation. With Gogny interaction we find that the equation 
of state (EOS) of asymmetric nuclear matter can be written as ^{p,5) w ^{p.O) + Esym{p)5'^ , in agreement with 
the empirical parabolic law found by all many body theories. Different calculations depending on the many body 
approaches and the interactions used, give very different behaviors for Egymip), especially at high densities. In 
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some, such as, the Relativistic Mean Field model [12] and the Brueckner-Hartrec-Fock [11], it is a continuously rising 
function of p. In others, such as, the Variational Many-Body Approach [13], it rises in the beginning and then begins 
to fall. Within the Hartree-Fock approach using all 86 Skyrme effective interactions widely used currently in nuclear 
structure studies, it was found that about 1/3(2/3) of them lead to symmetry energies in the first (second) category 
[14]. Gogny interaction with the default parameters has behaviour in the second category. It is interesting to see 
what causes the fall and how the parameters of the interaction can be changed minimally to alter this behaviour. In 
Fig. 2 we plot contributions to Esym{p) from the Gogny interaction and we show separately the term coming from (a) 
the direct term of the finite range part, (b) the exchange term of the finite range part and (c) the density dependent 
zero range part. It is the latter that causes the bending over. For example, if we choose the density dependent term 
to be S independent, Esym.{p) will continue to rise. 

Lastly, we write the potential energy density due to Gogny interaction in a form that is common practice in BUU 
literature. Thus 

V{p, 6) = -52.41(^) + 37.70(^)(52 + (1 - 6^) 



Po Po cr + 1 Po 

fr{f,P)fr'{fJ) [(-41.51 + 46.025,,,.) e-(^-^")'''?/4 



po ^ 

+ (-38.62 + 17.25(5,,/) e"'^-^')'"^/'*] dVV- (2.9) 

Here, as in previous work, all quoted numbers are in MeV, po = 0.16/m~'^. Also a = Q:(Gogny)+l=4/3. 

Given that the momentum dependence generated by the Gogny potential comes from two Gaussians, we proceed to 
find a simpler version for momentum dependent potential for asymmetric matter. We will extend the parametrisation 
of [3,4]. Subsequently we will find an even simpler version, the kind that was used in first applications of momentum 
dependent potentials for heavy ion collisions. 

III. A SIMPLE MOMENTUM DEPENDENT POTENTIAL FOR BUU CALCULATIONS 

The simplest generalisation of the potential energy density of Eq.(5.4) of [3] to asymmetric nuclear matter is 

v{p, s) = + ^p^s^ + 4t^(i - 

2po 2po <T + 1 p^ 



The parameter x is introduced to cover the largely uncertain behavior of nuclear symmetry energy Esym{p) as discussed 
in the previous section. For the choice of a; = 1 (same as in Gogny) in the term containig B, the symmetry energy 
will bend over beyond a density p; for the choice x — the symmetry energy will continue to rise with density. In the 
above, Ci/2,1/2 = (^-1/2,-1/2 = Qi/ce and Ci/2,-1/2 = ^^-1/2, 1/2 = Cuniike- In terms of interactions between like and 
unlike particles, the above equation is equivalent to 

v(p,5) = + ^{pI + pD + -§7^(1 - ^^') 

Po ^Po + 1 Po 

+ 1VC ,f f d'nd'n'MjPlMflll (3 2) 

PO ^^'^ J + ip-ff/A' ^^^^ 

where Ai = {Ay, + Ai)/2 and A2 = {Ai — Au)/2. The one-body potential needed for BUU computations is given by 
U{p,S,p,r) = AyP^ + aA + Bi^ra - xS') - x f /T 

Po Po Po £7 + 1 p5 dpr 



+ PO y ^ 1 + (p - p-) + PO y ^ 1 + {p- i?) ^ 



In the above r ^ r' and |^ = and |^ = 

' apn P OPp p 
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The constants appearing in Eqs.(3.1) and (3.2) will be fixed by ensuring that properties of cold nuclear matter are 
reproduced. There fri^jP) = 7^0(p/(t) — p)- The integration in Eq.(3.1) is facilitated by noting that for a fixed 
p= {pi — p2)/2, the centre of mass momentum can be integrated out to give 
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where qf = (pfir) +p/(t'))/2. 
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Fig. 3 Esym{p) ds a function of p for the two choices of x; 0(left panel) and 1 (right panel). 

For completeness, for cold matter, we write down the values of the integrals appearing in Eqs.(3.1) and (3.2) 



I fr{r,p)fr'{r.,p') 
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Similarly, the value of the integral in Eq.(3.3) is 
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Fig. 4 One body potential in cold matter for S ~ 0.4. As in fig. 3, here also the two panels are for x=0 and 1. 
Comparison with one-body potential as obtained with a Gogny interaction has also been shown. The solid line is the 
single particle potential for symmetric matter with MDI. 



We fix the force parameters by first optimally reproducing the variation of U{po, 5,p) with p with that obtained by 
using a Gogny force (Brueckner-Hartce-Fock calculation gives similar results [11]). This fixes A of Eq.(3.3). Other 
parameters are then fixed by values of saturation density (0.16 fm~^), binding energy (16 MeV) at saturation density, 
compressibility (« 210 MeV) of N=Z nuclear matter at saturation density and symmetry energy at TV = Z (« 30 
MeV). The values of the parameters for two choices (1) x=0 (denoted by MDI(O)) and (2) x=l (denoted by MDI(l)) 
in Eq.(3.2) as in Gogny (this causes the value of the symmetry energy to bend over as a function of density) and are 
given in the table. Fig. 3 shows the behaviour of symmetry energy for the two choices. In fig. 4 we show the one 
body potential in cold matter for 5 = 0.4. 

Before we close this section wc like to add that BUU calciilations with momentum dependent interactions arc only 
slightly more complicated in the asymmetric case compared to the case where no distinction is made between neutrons 
and protons. In numerical solutions of Boltzmann equations, phase space densities are simulated by test particles, 
characterised by a position and a momentum. Now there will also be a tag on their charge but no major addition to 
the codes [6] are needed. 



IV. REDUCTION TO THE GBD FORM 



Although we do recommend the full formalism of the above section be implemented, it is possible to reduce the 
above to a GBD (Gale, Bertsch and Das Gupta) form. A GBD potential, extended to asymmetric matter, already 

exists and is called BGBD [8]. There the extensions were made such that in the N = Z case one gets back exactly 
the original parameters [1]. The parameter x (eq.(3.1) was chosen to be 1/15. In contrast, here we have chosen 
the momentum dependence of the Gogny potential as a reference curve and chosen x at or 1. The values of the 

paramaters of GBD(O) {x=0), GBD(l) (.t=1) and BGBD (.t=1/15) arc given in the table. 
We write the potential energy density coming from the momentum dependent part as 
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The one-body potential generated by this piece of potential energy density for given p, 5, p and a given r is 



Umom{p,S,P,T) = 



+ 



Po 



h 



+ 



Pr 



Po 



P(f,p')dV 



1 + <P>r)VA', 
Pr' 



+ - 



(4.2) 



In the above equation, r' ^ r. As expected, the values of the constants in the force will have to be recalculated 
to reproduce the saturation properties. These are given in the table. The GBD potential U (p, S, p, r) is plotted in 
Fig. 5. This docs not track the Gogny potential as faithfully as the more sophisticated version of section III does. The 
reader might wonder why the value of A is significantly bigger in GBD as opposed to in MDI(O) and MDI(l). We 
have tried to fit the variation of U with k as obtained in Gogny potential (or the Brueckner-Hartree-Fock calculation) 
by adjusting the value of A. For p close pp, the contribution to U{p) in section III comes mainly from p' near pp 
whereas in GBD what counts is | p— < p >ave\ ■ since , <p >ave is zero, one requires a different value of A to mimic 
the variation with p. This point was not appreciated in [3]. 
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Fig. 5 One body potential in cold matter for 6 = 0.4. As in fig. 4: here also the two panels are for x=0 
1. Comparison with one-body potential as obtained with a Gogny interaction has also been shown. The solid line 
represents the single particle potential for symmetric matter with the GBD potential. 



V. MOMENTUM DEPENDENCE OF THE SYMMETRY POTENTIAL 



The single particle potentials derived in the previous sections can be used directly in transport model calculations. 
They combine the density, momentum and isospin dependences of both the isoscalar and symmetry potentials in a 
nontrivial way. In this section we evaluate the strength of the momentum dependence of the symmetry potential. To 
the leading order in 6, the single nucleon potential can be cast to the form 

Un/p{p,P,S) W Uo{p,p) ± Usym{p,P)S (5.1) 

in accordance with the Lane potential [17], where the ± sign is for neutrons and protons, respectively. Thus the 
symmetry potential can be evaluated from Usym{p,p) = iUn — Up)/ 26. Shown in Fig. 6 are the symmetry potentials 
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as a function of k for the three densities. It is seen that the Usyrn{p,p) is strongly momentum dependent for k < 5fmr^ 
in all models considered. Moreover, this dependence is particularly stronger at high densities. By construction, the 
results for Gogny and MDI(l) are very close. By comparing the results with a; = and a; = 1, it is seen that the 
momentum dependence is rather different mainly for p/po = 2. This is because the symmetry energies with a; = 
and X = 1 are significantly different only in the region of p > po as shown in Fig. 3. 
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Fig. 6 Dependence of (f/jv — Up)/25 on momentum for different potentials at different densities. All the figures are 
for 5=0.4 

To our best knowledge, the symmetry potential has been assumed to be momentum independent in all previous 
studies. Our results above indicate that this is only a good approximation at high momenta. At momenta less than 
about 1 GeV/c, the momentum dependence of the symmetry potential is important. 

VI. SUMMARY 

The goal of this work was to obtain parameters of a momentum dependent mean field potential which is applicable 
to highly isospin asymmetric nuclear matter but easy enough to use in a transport model calculation for heavy 
ion collisions at intermediate to high energies. We used the Gogny interaction as a guide. Published results of 
Brueckner-Hartree-Fock calculations were also used for choosing parameters. We have two versions. They both are 
quite flexible in the sense that parameters can be easily chosen according to that experimental data on binding energy, 
compressibility etc. In each of these versions we have proposed two sets of parameters: one, in which the symmetry 
energy continues to rise as a function of density and another one where the symmetry energy first rises and then falls 
off with density. Our opinion is that introducing momentum dependence in symmetry potential will not not make 
transport model simulations much harder or longer than they already are. Implementations of these potentials in 
BUU transport models are in progress. 

Transport models with the momentmn dependent symmetry potentials are more reliable tools for investigating the 
density dependence of nuclear symmetry energy and thus the EOS of neutron-rich matter. Knowledge on the symmetry 
energy Esym{p) is essential for understanding not only the structure of radioactive nuclei but also many key issues in 
astrophysics. For instance, the Esym{p) determines uniquely the proton fraction in neutron stars at /3 equilibrium [8]. 
A continuously rising symmetry energy leads to a growing proton fraction with increasing density, thus allowing for the 
fast cooling of protoneutron stars through the direct URCA process [15]. A falling symmetry energy at high densities 
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forbidcs the direct URCA process to happen; moreover, it favors the formation of pure neutron domains in the cores 
of neutron stars [16]. Nuclear reactions induced by neutron-rich nuclei provide a great opportunity to pin down the 
density dependence of nuclear symmetry energy. Among the experimental observables that have been found to be 
sensitive to the symmetry potential, the neutron/proton ratio of pre-equilibrium nucleon emissions, neutron-proton 
differential flow and correlation functions, as well as the proton elliptic flow at high transverse momenta are expected 
to be most sensitive to the momentum dependence of the symmetry potential. These observables will be studied 
with the improved BUU transport models using the momentum dependent symmetry potentials. These results will 
be reported in a forthcoming publication. 
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TABLE I. The values of the parameters of different interactions. Also gives the values K in the respective parametrisations. 
The saturation density for all parametrisation is 0.16/m~^, except for BGBD(1/15) it is 0.163/m~^. The binding energy and 
the total symmetry energy including contribution from kinetic part are -16 MeV/A and 31.623 MeV/A respectively, in each 
case. 



Parameters 


MDI(O) 


MDI(l) 


GBD(O) 


GBD(l) 


BGBD(1/15) 




-95.98 


-187.27 


-109.85 


-299.69 


-192.0 




-120.57 


-29.28 


-191.30 


-1.46 


-96.0 


B 


106.35 


106.35 


205.66 


205.66 


203.3 


Cu 


-103.40 


-103.40 


-118.80 


-118.80 


-84.53 


Ci 


-11.70 


-11.70 


-26.26 


-26.26 


-65.472 


a 


4/3 


4/3 


7/6 


7/6 


7/6 


A 


1.0p/°) 


1.0p/°) 


1.5p/°' 


1.5p/(») 


1.5p/°' 


K 


210.68 


210.68 


214.68 


214.68 


215.0 
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